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The  background  of  n.ierostructure  and  its  effect  on  sound  propagation 
is  reviewed.  The  expendable  bathythermograph  (XBT)  system  and  current 
digitizing  procedures  are  evaluated  to  determine  their  adequacy  in 
investigating  vertical  thermal  microstructure .  Several  methods  of 
microthermal  analysis  are  proposed.  Temperature  gradient  versus  depth 
plots  are  used  to  analyze  data  from  several  ocean  regions  -  Atlantic 
Coastal,  Andaman  Sea,  Eastern  North  Pacific,  and  the  northern  boundary 
of  the  Gulf  Stream.  Analysis  of  these  data  suggest  a  water  mass  depend¬ 
ence  of  the  thermal  microstructure.  The  significance  of  the  small-scale 
structure  to  Naval  operations  is  discussed. 
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I.  SIGNIFICANCE  OF  MICROTHERMAL  STRUCTURE 


Nava]  interest  in  the  thermal  structure  of  the  ocean  is  presently 
centered  on  gross  features  such  as  sea  surface  temperature,  mixed  layer 
depth,  thermocline  gradient,  and  depth  to  the  bottom  of  the  thermocline. 
Some  attention  has  been  given,  especially  in  the  Canadian  Navy,  to  the 
transients  in  the  mixed  layer  [Tully,  1964]  and  their  possible  associ¬ 
ation  with  internal  waves.  The  above  gross  features  are  used  primarily 
as  a  means  of  classifying  thermal  structure. 

Technological  developments  have  led  to  measurements  which  indicate 
the  thermal  structure  is  made-up  of  layers  of  varying  thickness;  from 
less  than  a  centimeter  to  tens  of  meters  [Stommel  and  Federov,  1967; 
Osborn,  1969;  for  example]. 

Of  the  variables  affecting  sound  velocity  in  the  ocean,  temperature 
has  the  greatest  effect.  The  presence  of  variations  in  temperature,  or 
temperature  gradient,  causes  sound  velocity  fluctuations  to  occur. 
Research  on  the  effects  of  the  thermal  layers  on  sound  in  the  ocean  show 
that  further  investigation  in  theory  and  experiment  is  required.  Small 
scale  sound  velocity  variations  cause  irregular  bending  of  sound  rays. 
The  microstructure  can  affect  attenuation,  forward  and  backward  scatter¬ 
ing,  reverberation,  reflection,  and  refraction  which  can  converge  or 
diverge  the  energy  in  a  sound  beam  [Meyer  and  Romberg,  1963].  The 
physical  features  of  the  microthermal  structure  must  be  adequately 
described  before  their  effects  upon  acoustic  propagation  can  be  fully 
understood . 
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The  purpose  o f  this  thesis  is  to  review  the  present  knowledge  con¬ 
cerning  microstructure  as  measured  by  expendable  bathythermographs  (XBT) , 
to  evaluate  the  capability  of  the  XBT  in  measuring  microthermal  structure, 
to  develop  a  scheme  of  quantitatively  examining  the  microthermal  structure, 
and  to  teat  this  scheme  on  real  data. 
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II.  REVIEW  OF  PRESENT  KNOWLEDGE 


A.  MXCEOTHERMAL  STRUCTURE  IN  THE  OCEAN 

Prior  to  1937,  vertical  theraal  structure  was  measured  at  discrete 
points  using  reversing  thermometers.  This  method  provided  a  point-to- 
point  analysis  of  the  vertical  structure  at  one  location.  Even  using 
point-to-point  measurements ,  it  was  clear  that  the  temperature  vari¬ 
ability  decreased  with  depth  in  most  regions.  In  fact,  the  most  common 
analysis  of  the  temperature  and  salinity  data,  T-S  correlation,  usually 
ignored  the  upper  100  meters  or  so.  However,  this  same  analysis  showed 
that  in  many  ocean  regions,  in  deep  water,  (e.g.,  below  1000  to  2000  m), 
the  variability  was  on  the  order  of  the  instrument  error  [Sverdrup, 

*  Johnson  and  Flamming,  1942].  This  result  was  not  unexpected  because 

\  the  gradients  in  the  deep  ocean  became  smaller.  The  discrete  measure¬ 

ments  of  the  Nansen  cast  are  not  suitable  for  studying  small  scale 
variability. 

The  development  of  the  mechanical  bathythermograph  (MET)  by  Spilhaus 
[1938]  provided  the  first  convenient  means  of  measuring  a  continuous 
profile  of  temperature.  The  mechanical  bathythermograph  was  widely  used 
during  World  War  II.  A  large  number  of  be thy thermograms  with  double 
traces,  which  cou*.*'  not  be  explained  as  instrument  errors,  were  dis¬ 
covered  [Urick  and  Searfoss,  1948].  The  double  traces  indicated  that 
small-scale  horisontal  and/or  vertical  temperature  variations  occurred 

? 

between  the  lowering  and  raising  of  the  MBT.  Urick  and  Searfoss  [1949] 
recognised  that  this  small-scale  microstructure  could  be  significant  to 
sound  propagation,  depending  on  the  site  and  shape  of  the  inhomogeneity. 
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Horizontal  mlorothfcrmai  structure  was  measured  by  Bolder  in  1944 


using  a  thermopile  mounted  on  a  submarine.  He  detected  horizontal 
temperature  differences  of  0.C2?  over  distances  of  about  10  m  (Smith, 
196?}.  Measurements  of  these  horizontal  temperature  variations  were 
also  istde  by  Urick  and  Searfoss  (1948,  1949}  in  July  of  1948  near  Kay 
Uhat,  Florida.  With  a  thermocouple  mounted  on  the  conning  tower  and, 
later,  its  bow,  a  submarine  proceeded  at  a  constant  spaed  and  dapth, 
while  tha  probe  measured  the  water  temperature  to  a  resolution  of 
O.OdSC. 


0EPTH-  25 ft 


Figure  l.  Hicrothermal  Variations  at  Various  Depths  , 

[Urick  and  Searfoss,  1948). 

Urick  and  Searfoss  (1949)  calculated  the  peak-to-peak  root-mean- 
squarc  (rms)  variation  for  peaks  in  excess  of  20  percent  of  the  total 
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variation  in  a  record  calling  this  the  root -mean-square  maximum  (rmsm). 

A  significant  difference  between  the  value  of  rmsm  in  the  mixed  layer 
and  in  the  thermocline  was  observed.  In  the  mixed  layer  a  mean  rmsm 
value  of  0.013C  was  observed.  In  the  thermocline,  the  mean  value  was 
0.450C.  The  rmsm  increased  with  depth  in  the  mixed  layer  and  decreased 
with  depth  in  the  thermocline.  Applying  Taylor's  hypothesis,  the  mean 
horizontal  thermal  "patch  size"  was  found  to  vary  from  5  to  30  m  in  the 
mixed  layer  and  from  30  to  200  m  in  the  thermocline  region  [Urick  and 
Searfoss,  1949]. 

Liebermann  [1951]  mounted  a  sensitive  platinum-resistance  thermom¬ 
eter  and  a  thermocouple  on  the  periscope  of  a  submarine.  The  platinum- 
resistance  element  could  resolve  spatial  temperature  variations  on  the 
order  of  10  cm  in  horizontal  distances  as  the  submarine  proceeded  at  a 
constant  speed  of  a  few  meters  per  second  at  a  fixed  depth  (Figure  2). 

To  characterize  the  temperature  fluctuations,  he  used  the  autocorrelation 
function.  "Patch  size"  was  defined  as  the  distance  of  the  microthermal 
structure,  within  which  the  temperature  variations  maintained  coherence. 

The  function  is  given  by: 

00 

R  (o)  -  f  [T  (x)-T  (xfp0)]  dx.  (1) 

0 


His  results  could  be  fit  equally  well  by  an  exponential, 

a  Co)  -  e-°/a 

or  the  Gaussian, 

R  (o)  -  e‘(p/a) 

where  ? 


(2) 


(3) 
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T  (x)  *  tempcratuta  at  soma  point  x; 

T  (x+Pq)  “  temperature  at  another  point  displaced  a  distance  p^ 
from  x ; 

o  *  autocorrelation  distance; 

<p>  *  a  *  "patch  site". 

Whereas,  the  exponential  function  fitted  the  data  points  and  decayed 
rapidly  as  p  increased  (Pigure  2),  it  did  not  have  a  aero  slope  at  the 
origin,  thus,  the  exponential  function  was  not  appropriate  for  very 
small  correlation  distance  (e.g.,  a  few  cm  or  less)  because  of  heat 
conductivity  in  sea  water.  He  proposed  the  Gaussian  function  to 
correct  this  difficulty,  reasoning  that  thermal  conductivity  would 
prevent  any  discontinuity  in  the  temperature  field  (Liebermann,  1951]. 

Liebermann  assumed  that  the  thermal  microstructure  represented  a 
cloud  of  spheres  with  a  Rayleigh  distribution  having  a  mean  else  value 
of  the  correlation  distance  and  therefore  yielding 

R  (p)  ■  R  (a)  ■a"1.  (4) 

Using  this  definition,  the  mean  slae  of  the  inhomogeneities,  a,  was 
found  to  equal  60  cm  in  the  ocean  (Figure  2). 

Piip  [1961]  investigated  the  sound  velocity  microstructure  in  the 
sound  channel  off  Bermuda  with  a  NBS  velocimeter.  It  was  noted  that  a 
number  of  thin  steps  remained  over  periods  of  hours,  some  oscillating 
30  m  in  depth  in  a  three-hour  period.  At  the  same  station  and  same 
depth,  sound  velocity  changes  of  one  m/sec  in  a  few  hours  were  common 
and  "nearly  the  rule".  The  sound  velocity  axis  was  observed  to  change 
150  m  in  depth  and  3.5  m/sec  over  a  three-month  period.  Pipp  showed 
that  the  base  of  the  main  thermocline  (about  1000  to  1300  m)  exhibited 
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relatively  large  and  rapid  change*  in  sound  velocity  in  an  inhomogeneous 
region  in  a  constant  state  of  fluctuation. 


•  meters 


\ 


R(p) 


p,  cm 


Figure  2.  Temperature  Fluctuations  and  Associated  Correlation  Function 

(After  Llebermann) 

8mall-scal*  temperature  and  salinity  structure  was  measured  with  a 
salinity i  temperature,  and  depth  recorder  (STD)  by  Stomnel  and  Federov 
[1967]  near  Timor  and  Mlndinao  during  July  of  1965.  Their  measurements 
revealed  that  the  vertical  microstructure  actually  consisted  of  very 
thin  high-gradisnt  layers.  The  thin  layers  were  separated  by  thicker 
layer-  of  nearly  homogeneous  (isothermal  or  isohailne)  water,  i.e., 
step- like  in  appearance  on  their  STD  traces.  By  following  a  discrete 
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feature  from  one  cast  to  the  next,  homogeneous  layers  were  found  to  be 
continuous  from  two  to  20  km  horiz9ot*lty  and  to  range  in  thickness  from 
two  to  40  m  (Figure  3).  A  boot-like  feature  was  found  at  the  bottom  of 
the  mixed  layer  (Figure  4).  At  a  depth  of  200  to  300  m,  the  mean  thermal 
gradient  was  about  0.040C/m,  J.n  the  thin  steps,  the  gradient  was  less 
than  0.003C/m,  and  in  the  homogeneous  steps,  the  thermal  gradient  was 
greater  than  0.350C/m.  Stommel  and  Federov  concluded  that  the  ocean 
consisted  of  nearly  homogeneous  layers  separated  by  steps  with  extremely 
high  thermal  gradients. 

Tate  and  Howe  [1968],  in  the  summer  of  1966  in  the  Northwestern 
Atlantic,  used  an  STD  to  look  at  the  thennocline  stratification  in  the 
Mediterranean  Water  intrusion  at  a  depth  of  1280  to  1500  m.  They  found 
layering  with  steps  varying  on  the  order  of  0.17  to  0.35C  in  tempera¬ 
ture,  from  0.02  to  0.055  o/oo  in  salinity,  and  from  thickness  of  15  to 
30  m  (Figure  5).  A  step's  average  temperature  change  of  0.25C  produced 
a  3  change  of  -0.044,  while  the  average  salinity  step  of  0.044  o/oo 
produced  a  ot  change  of  +0.039.  The  resulting  change  of  o  (-0.005) 
gave  a  slightly  stable  layer  system  in  a  situation  which  was  clearly 
one  of  opposing  salinity  and  temperature  density  gradients. 

Cooper  and  Stommel  [1968],  using  an  STD,  found  the  main  thermocline 
in  the  Sargasso  Sea  to  consist  of  regular  salinity  and  temperature  steps 
of  homogeneous  layers  three  to  five  meters  thick,  separated  by  trans¬ 
itional  layers  10  to  15  meters  thick,  with  temperature  changes  of  0.30 
to  0.5QC  and  salinity  changes  of  0.04  to  0.10  o/oo  (Figure  6).  The 
horizontal  extent  of  these  steps  was  between  400  to  1000  meters.  The 
main  thermocline  contained  up  to  100  steps  on  each  cast. 
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a. 


b. 


Figure  3. a.  Photographic  copy  of  actual  STD  trace  (not  retouched)  of 
station  21  off  Hlndinao.  Depth  scale  is  500  m,  each  small 
division  being  5  m.  The  salinity  curve ,  labeled  S,  has 
the  smallest  division  being  0.05  o/oo.  The  temperature, 
labeled  T,  has  a  smallest  division  of  0.05C.  The  lamina 
B  is  so  labeled.  Because  the  pens  cannot  both  traverse 
the  same  path  mechanically,  the  temperature  curve  is 
offset  downward  by  5  m  and  always  reads  5  m  too  deep. 

The  salinity  scale  is  not  offset  in  depth.  Uniform 
corrections  for  salinity  of  -0.9  o/oo  and  temperature  of 
0.07C  have  not  been  applied. 

Figure  3.b.  A  composite  chart-graph  for  lamina  3  off  Mindinao.  Each 
square  is  0.05C  wide  (8.5  to  9.0C)  and  50  m  high.  The 
small  arrow  indicates  p  depth  of  400  m  in  each  square. 
Numbers  in  the  lower  le^'t  hand  corners  are  the  station 
numbers.  The  distances  from  stations  18  to  25  and  26  to 
35  are  about  10  nm.  For  lamina  B,  the  large  geographical 
extent  of  this  small  lamina  is  self-evident. 

[Stommel  and  Federov,  1967] 
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Figure  4.  From  the  soundings  in  the  mixed  layer  off  Timor,  gross  dif¬ 
ferences  in  the  properties  in  different  locations  are  self- 
evident.  There  is  considerable  variation  in  vertical 
structure  but  each  sounding  shows  a  nearly  homogeneous 
region  near  the  surface  (station  13)  but  some  stations 
(like  7)  have  shallow  Irregular  layers  overlying  the  mixed 
leyer. 


[Stoomel  and  Federov,  1967] 

T&S 


Figure  5  A  ^vplral  step  in  Mediterranean  Water  (Talt  and  Howe,  1968). 
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Woods  [1968,  1968a]  investigated  the  micros true ture  in  the  summer 
thermocline  off  Malta  using  a  temperature-gradient  meter  for  obtaining 
the  thermal  gradient  between  two  thermistors  set  50,  25  or  10  cm  apart 
on  a  vertical  staff.  The  resolution  of  the  thermal  gradient  was  within 
0.01C  per  thermistor  separation  distance.  An  additional  thermistor 
recorded  temperature.  Woods  plotted  profiles  of  the  thermal  gradient 
and  temperature  versus  depth  (Figures  7,  8  and  9).  Using  time-lapse 
photography ,  Woods  obtained  pictures  of  streaks  left  in  the  wakes  of 
free-dropping  dye  pellets  for  shear  measurements.  A  dye-packet  array 
tied  between  moored,  submerged  floats  injecting  dye  into  selected  levels 
of  the  thermocline,  indicated  that  the  thermocline  was  divided  into 
layers  a  few  meters  thick.  These  layers  were  characterised  by  weak 
thermal  gradients,  less  than  O.OOlC/m,  and  by  turbulence  with  a  mean 
velocity  of  less  than  one  mm/sec.  The  layers  were  separated  by  sheets 
only  centimeters  thick.  (Sometimes,  he  detected  sheets  up  to  a  meter 
thick  but  they  appeared  as  aggregates  of  the  10  cm  thick  sheets.)  The 
sheets  were  characterized  by  strong  thermal  gradients  up  to  0.05C/cm 
and  little  or  no  turbulence.  Mature  sheets  contained  thin  laminar-flow 
zones  in  their  centers.  Neighboring  layers  had  dissimilar  properties 
indicating  different  salinities,  velocities,  and  turbidities.  This 
suggested  different  trajectories  and  a  very  slow  exchange  of  properties 
across  the  intervening  sheets.  If  turbulence  was  small,  flux  exchange 
was  small;  therefore,  the  exchange  had  to  rely  on  molecular  processes, 
which  were  very  slow  indeed. 

Individual  layers  were  followed  horizontally  for  distances  up  to 
72  km  by  soundings  from  high-speed  launches.  Layers  were  identified  on 
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Figure  8 


Temperature  and  Temperature-Gradient  Soundings  for  Maltese 
Waters,  14  September  1967  [Woods,  1968a]. 

Temperature  0918-0919  local  time,  (i.e.,  GMT  +  1). 

Temperature  gradient  1132-1145  local  time. 
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successive  soundings  for  many  hours,  thus  confirming  the  results  of 
Stomiei  end  Federov. 

The  dye  experiments  showed  thet  in  a  thermocline  sheet  there  was  a 
central  laminar-flow  region  only  a  few  centimeters  thick.  Dye  was  spread 
in  this  region  by  internal  waves  to  form  a  "broad  carpet  of  color".  This 
colored  sheet  made  it  possible  to  photograph  the  internal  waves  as  they 
crossed  it.  Wave-lengths  between  5  and  10  m  were  dominant  and  appeared 
in  long,  coherent  trains,  each  associated  with  a  single  sheet.  They  had 
heights  up  to  1,0  o,  a  phase  speed  of  2.5  cm/sec  and  periods  of  about 
5  min.  These  waves  only  affected  a  single  sheet,  being  much  shorter  than 
the  long  internal  waves  which  move  the  entire  thermocline. 

A  sheet  was  defined  as  stable  when  it  appeared  smooth  everywhere ; 
this  indicated  laminar  flow.  When  a  wave  passed  over  a  sheet,  shear  had 
its  greatest  effects  at  the  crest  and  trough.  Wave  shear  and  drift 
shear  interacted  either  cancelling  or  reinforcing  each  other  causing  a 
sheet  to  become  unstable  and  patchy. 

Wavelets,  about  75  cm  in  length,  caused  banding  to  appear  parallel 
to  the  crests  and  troughs.  These  wavelets  generally  formed,  grew,  and 
broke  in  a  period  of  less  than  two  minutes.  Upon  breaking  in  classical 
fashion,  a  second  smaller  breaker  was  generally  thrown  forward.  The 
wavelets  had  a  maximum  height  of  about  20  cm.  After  breaking,  a  turbu¬ 
lence  region  remained  for  about  5  minutes,  and  the  scar  in  the  sheet 
remained  for  several  hours.  These  scars  were  a  common  feature  in 
thermocline  sheets.  The  thinnest  sheet  observed  to  become  unstable  had 
a  thickness  of  3  cmj  the  thickest  was  30  cm.  The  great  majority  were 
between  8  and  15  cm. 
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Transient  thermoc i tn.ss  [Tully,  1964]  which  occur  sporadically  between 
the  surface  and  the  tup  of  the  thermoc line  were  considered  aa  further 
examples  of  sheets.  Woods  divides  the  thermae  line  into  a  half-dosed  or 
so  layers  of  low  shear  and  moderate  thermal  gradient  separated  by  thin 
llminar  flow  sheets  on  intense  shear  and  thermal  gradient.  He  discussed 
the  important  aspect  of  heat  transfer  through  such  a  layered  ocean.  He 
suggested  that  the  molecular  heat  transfer  through  the  high  gradient 
stable  layers  was  negligible,  but  that  heat  transfer  does  take  place 
through  the  turbulent  scars  from  one  layer  to  another  as  an  effect  of  the 
breaking  of  the  internal  waves. 

The  use  of  the  Ramsey  SVTP  (sound  velocity,  temperature,  and  pressura) 
probe  provides  a  digital  output  of  these  parameters.  This  instrument  has 
an  accuracy  of  4;  0.01C  and  +  0.25  percent  of  depth  with  a  precision  of 
+  0.08  m/sec  in  sound  velocity  [Lovett,  1968].  In  view  of  the  plus-or- 
minus-one  count  for  electronic  counters,  small*scale  features,  i.e., 
less  than  4;  0.01C  and  +  1  si  in  depth,  are  ambiguous.  Lovett  (1968] 
discussed  the  thermal  gradients  determined  with  this  instrument  (Figure 
10).  The  10  m  depth  increments  he  used  eliminated  some  microstructrre ; 
however,  using  5  m  increments  introduced  random  readout  errors, 

Inghtm  [1968]  used  an  STD  probe  to  investigate  the  mixed  layer  in 
tropical  waters.  The  vertical  salinity  and  temperature  structures 
forming  the  "Isothermal"  and  "isohaline"  layers  were  generally  analogous 
in  the  near-surface  layers.  The  difference  in  the  depths  of  these  two 
layers  exceeded  7.5  m  in  45  percent  of  the  casts  and  17.5  m  in  31  percent 
as  sh^wn  by  Figure  il  [Ingham,  1968].  Microstructure  effects  on  sound 
velocity  in  the  near-surface  layers  would  depend  on  salinity  in  some  ocean 
regions ., 
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figure  10.  Temperature  Gradient  per  10  a,  Keen  Gradient  per  100  *,  end 
Average  Deviation  from  Keen  versus  Depth  team  tommy  Probe 
Date  [Lovett,  1908]. 


figure  11.  Histogram  of  Depth  Differences  AZ  of  the  Isothermal  and 
laohaline  layer  Depths  [Ingham,  1968]. 
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Pingree  [1969]  using  an  STD,  recorded  analogue  signals  on  tape.  He 
made  fine-scale  plots  of  salinity  and  temperature.  In  particular,  he 
Investigated  the  stability  of  the  deep  water,  i.e.,  about  500  to  2000  m 
(Figure  12), 

Osborn  [1969]  using  a  free-falling  temperature  gradient  meter 
obtained  excellent  profiles  of  temperature  gradient  (Figure  13).  Off 
San  Diego,  the  horizontal  extent  of  the  layers  was  greater  than  750 
meters  in  the  seasonal  thermocline,  but  only  a  few  hundred  meters  at 
depths  greater  than  400  meters.  He  noted  the  layers  were,  in  fact, 
made  up  of  even  finer  layers  with  a  thickness  of  only  a  few  centimeters. 
His  average  sheet  was  about  one  meter  thick. 

Recent  work  by  Neshbya,  Neal,  and  Denn^r  [1969]  using  the  expendable 
bathythermograph  manufactured  by  General  Motors  Defense  Laboratories, 
obtained  records  of  temperature  steps  (Figures  14,  15  and  16).  While 
on  Ice  Station  T-3  in  the  Arctic  (near  84-38N,  128-22W) ,  they  lowered 
XBT  probes  by  hand.  The  Ice  Station  results  are  the  first  XBT  micro¬ 
structure  records  without  ship  motion. 

B.  APPLICATIONS 

Urick  and  Searfoss  [1949]  calculated  theoretically  that  a  sound  ray 
would  be  bent  through  an  angle  of  only  3.7  sec  by  the  average  inhomo¬ 
geneity  in  the  mixed  layer  but  as  much  as  2.1  min  of  arc  in  the 
thermnc line . 

Sheehy  [1950]  measured  scnar  signal  fluctuation  as  functions  of  depth 
and  range.  He  found  the  coefficient  of  variation,  defined  as  the  stand¬ 
ard  deviation  of  the  amplitudes  expressed  as  a  percentage  of  the  mean 
amplitude,  varied  as  t he  square  root  of  range  and  did  not  appear  to  be  a 
function  of  depth. 
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Figure  12 .  T-S  Curve  f 
Circles  Rap 
1969], 
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150m 


(a)  Typical  Temperature 
Profile  Section 


IfxMJ 


figure  14.  Vertical  Profile  of  Temperature  Under  Ice  Station  T-3 
(•4-3811,  128-22W)  19  March  1969.  [Neahbya,  Neal  and 
Danner,  1969]. 
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Figure  16.  tlM  krtu  with  High  Anpllficatloo  Taken  fro*  T-3  in  March 
1969.  Tempera  tura  lacraaaea  fro*  right  to  left.  [Raahbya , 
Maal  and  Daanar,  1969]. 
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Liebersutnn  [1951]  described  the  influence  of  the  micros true Cure  on 
acoustic  propagation  Be  wrote  the  sharpness  of  temperature  boundaries 
determines  influence  of  the  microstructure  on  acoustic  propagation. 

When  temperature  changes  occur  within  a  distance  of  less  than  one 
acousti  wavelength,  the  temperature  boundaries  may  reflect  sound 
energy  depending  upon  the  sise  and  shape  of  the  "patch"  The  tempera* 

i 

ture  microstructure  may  reflect,  refract,  scatter,  or  focus  acoustic 
waves,  The  microtherma 1  structure  was  considered  to  be  made  up  of 
randomly  slsed  "spheres",  spaced  randomly  throughout  tha  medium,  with 
e  scattering  cross  section  dependent  upon  the  autocorrelation  of  tha 
temperature  inhomogeneities,  R  (o) .  Comparing  derived  scattering  cross 
section  with  actual  scattering  in  the  ocean  (i.e.,  volume  reverberation), 
Liabermann  concluded  the  observed  reverberation  was  higher  than  one 
would  expact  from  thermal  inhomogeneities  He  attributed  volume 
reverberation  to  biological  scatters. 

Liabermann  found  that  refraction  affects  of  tha  inhomoganeltles 
resulted  in  warping  of  tha  wavs  fronts,  and  caused  tha  wave  front  to 
be  a  complex  surface  exhibiting  both  convergence  and  divergence.  Tha 
temporal  motion  of  the  inhomogene  it ia a  resulted  in  sound  "scintillation" 
at  tha  receiver 

Mlntrer  [1953]  found  Sheehy'a  coafficiant  of  variation  to  be  depend* 
ant  upon  the  three-halves  power  of  the  range  rather  than  the  one*half 
power  of  range . 

Skudrzyk  [1957]  discussed  the  focusing  and  defocuaing  affect  within 

2 

the  focusing  range  rQ  given  by  rQ  •  ka  where  a  la  tha  "patch"  radius 
and  k  is  the  acoustic  wavenumber  At  ranges  greater  than  r^,  in  the 
Interference  range,  phase  interference  results  (Figures  17  and  18). 


Figure  17.  Illustration  of  the  Various  Type*  of  Acoustic  Scattering 
Oua  to  a  Thernal  Fateh  in  tha  Sea  (Whitemarsh,  Skudrsyk 
and  Urick,  1957] . 


Figure  18,  Tha  Riaaa  Angle  Between  tha  Incident  and  tha  Scattered 
Praeeure  and  tha  Fhyaical  Nature  of  tha  Scattering  a?  e 
Function  of  the  Sanga  [Skudrsyk,  1957], 
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Using  this  concept,  Skudrsyk  (Figure  19)  show*  th*  focusing  range  r^  as 
a  function  of  frequency  and  patch  else.  (For  example,  a  3G  MU  frequency 
sonar  has  a  focusing  range  of  75  meters,  if  the  patch  slse  is  assuamd  to 
be  60  centimeters.)  For  finding  patch  slse,  Skrudxyk  used  the  Urlck- 
Sear f oss  (1948)  data  to  show  patch  else  was  twica  the  depth  (Figure  20). 
Further,  Skudrsyk  started  theoretical  investigation  with  rpheroidal 
patches,  which  he  found  to  have  a  considerably  stronger  focusing  effect 
than  the  spherical  patches. 

Whlteoarsh,  Skudrsyk,  and  Urlck  (1957)  did  additional  work  with  the 
spheriodal  theory.  The  correlation  functions  previously  used  were  dis¬ 
cussed.  It  was  concluded  that  the  present  correlation  functions  ware 
not  suitable  and  either  had  to  be  replaced  or  redefined.  The  author* 
argued  if  the  temperature  microstructure  was  produced  by  turbulent 
motion  in  the  sea,  then  the  Kolmogorov  Law  for  atatistical  aquiilbrlum 
must  hold.  This,  thay  show,  was  in  axcallant  agrt ament  with  thair  data. 
Landing  support  to  continued  studies  such  as  are  presented  in  thia 
thesis,  their  final  conclusion  states  that  the 

"...  greatest  gap  stems  to  lie  in  the  nlcrothamei 
statistics  of  the  ocean.  If  these  wars  known  batter 
than  they  are  now,  apparently  the  amplitude  statistics 
of  the  fluctuations  of  sound  transsiltted  through  the  . 
set  could  be  predicted  roughly,  but  with  some  assurance." 

Works  by  Peterson  (1965],  Barakus  (1968]  and  DINapoll  (1969)  have 
considered  the  normal  mode  approach  to  the  interference  range  solution 
for  the  temperature  inhomogeneities  and  have  also  investigated  the 
effects  of  internal  waves  on  sound  in  stratified  mediums. 


Hfhitemarsh,  D.  C.,  Skudrsyk,  E.,  and  Urick,  R.  J. ,  "Forward  Scattering 
of  Sound  in  the  Sea  and  Its  Correlation  with  the  Temperature  Micro¬ 
structure,"  J,  Acous .  Soc .  Am. .  v.  29,  p.  1*41,  October  1957. 
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Figure  19.  Focusing  Xenge  (meters)  a*  •  Function  of  Frequency  and 
Patch  Site  (meters)  [Skudrsky,  1937]. 
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DEPTH  (ft) 


Flgu  re'  20,  Fateh  Diameter  Veraua  Dtpth  (after  Skudrtyk*  1957], 


C.  INADEQUACIES  OF  THE  PRESENT  THEORY 

Since  the  thermal  microstructure  is  not  related  simply  to  the  gross 
thermal  features  measured  on  MBTs,  there  is  presently  no  easy  way  to 
predict,  on  a  day-to-day  basis,  the  effects  of  microthermai  structure 
on  the  propagation  of  sound  in  the  sea.  The  theories,  so  far,  have 
developed  in  two  classes. 

The  first  theories  (Liebermann,  1951;  Mintser,  1953,  1953a,  1954; 
and  Skudrsyk,  1957]  concerned  spherical  and,  more  recently,  spheriodal 
shaped  inhomcgcneitles  throughout  the  ocean.  Part  A  above  has  dis¬ 
cussed  the  present  knowledge  of  microstructure  showing  that  the  long 
thin  sheets  and  layers  do  exist  as  microthermai  structure.  These  first 
acoustical  theories  do  not  apply  to  what  is  now  described  as  the  micro- 
structure  of  the  ocean. 

The  theoretical  treatment  by  Whitemarsh,  Skudrsyk  and  Urick  [1957], 
considering  sperlodal  patches  and  Kolomogorov's  turbulence  theories, 
could  not  find  any  real  theoretical  explanation  of  what  happens  in  the 
sea  when  sound  passes  through  an  Inhomogeneity.  The  investigations 
using  normal  modes  have  not  produced  any  practical  operational  results. 

Lee  [1961]  using  ray  path  theory,  demonstrated  the  effect  of  a 
typical  slnusiodal  internal  wave  o  sound  intensity  in  a  three-layered 
model.  For  the  medium  with  an  internal  wave,  intensity  variations  up 
to  22  db  were  found  in  distances  of  less  than  one  Internal  wavelength, 
as  shown  by  Figures  21  and  22,  In  contrast,  the  same  medium,  without 
an  Internal  wave,  produced  a  5  db  variation,  including  a  2.5  db 
spreading  loss,  over  the  same  distance. 

Another  study  investigation  the  effects  of  sound  velocity  micro¬ 
structure  was  conducted  by  FNWC  [Wolff,  Tatro  and  Megehee ,  1967],  A 
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figure  21.  Diagreu  of  Sound  lay#  through  •  Modlua  Which  ho*  on 
Internal  Wave  or  tho  Theraocllna  {loo,  1961]. 


mumo  intimity  nno  hoi 


ct~> »*»  racsai»»o  ezzzao**?.* 


Figure  22.  Sound  Level  in  tho  Hediua  with  on  Internal  Wove  on  the 

Thermocllne.  Tho  db  reference  level  la  that  corresponding 
to  e  aound  level  of  60  db  et  1  toot  frem  the  directional 
•ource  (X  equals  0,  Z  equals  10  ft)  along  tho  horliontal. 
The  field  la  contoured  at  Intervale  of  2.5  db  [Lee,  1961]. 


sophisticated  ray-tracing  program  [Ayers,  Wolff,  Carstensen  and  Ayers,  1966] 
was  usad  for  tha  analysis.  Three  comp lata  ray  traca  runs  war a  made  to 
tha  first  convsrganea  sons  using  120,  54,  and  20  dlgltlsad  sound  velo¬ 
city  points  from  one  SVTF  traca  containing  sound  velocity  micros true tura . 

Tha  120-polnt  sound  velocity  profile  preserved  all  tha  "wiggles"  or 
microstructure,  tha  54-point  profile  usad  values  at  avary  50  mater  of 
depth,  and  the  20-polnt  profile  usad  tha  standard  oceanographic  depth 
values  for  sound  velocity  (Figure  23).  Tha  range  to  tha  first  conver¬ 
gence  aone  agreed  within  150  yards  for  all  three  cases.  Based  on  the 
ray  trace  runs  only,  it  was  concluded  the  sound  velocity  microstructure 
was  "essentially  self -cancel ling"  and  did  not  affect  the  convergence 
sone  parameters  [Wolff,  Tatro,  and  Meghee,  1967].  Propagation  loss  and 
reverberation  effects  were  not  reported. 


Figure  23.  A  Portion  rf  the  Three  Dlff«r*'ut  Sound  Velocity  Profiles 
Used  by  FNWC  for  Ray  Tracing  [After  Wolff,  et  *1.1. 
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III.  MEASURING  M1CR0THBRMAL  STRUCTURE 


A.  THE  BATHYTHERMOGRAPH 

1.  Mechanical  Bathythermograph 

At  the  present  time  temperature  is  the  most  readily  observed 
variable  in  the  upper  few  hundred  meters  of  the  ocean.  A  major  con* 
tributor  to  this  situation  is  a  singularly  successful  oceanographic 
instrument  known  as  the  mechanical  bathythermograph  or  MBT.  Subsequent 
to  its  introduction  by  Rossby  and  Montgomery  [1934]  and  perfection  by 
Spilhaus  [1938],  the  MBT  has  become  one  of  the  basic  instruments  of 
almost  all  oceanographic  agencies  because  of  its  simplicity  of 
operation,  reliability  and  relatively  low  cost.  The  result  has  been 
that  a  large  number  of  bathythermograms  have  been  collected. 

Researchers  during  World  War  II,  using  the  MBT,  observed  micro* 
thermal  structure  on  their  traces.  No  theory  had  been  formulated  which 
took  into  account  the  effects  caused  by  mlcrothermel  structure  on  under¬ 
water  sound  transmission.  It  was  recognised,  however,  that  temperature 
mlcroatructure  must  be  considered  in  any  explanation  of  observed  under¬ 
water  sound  transmission  [National  Defense  Research  Committee,  1946]. 

2.  The  Expendable  Bathythermograph 
a.  Capabilities 

Presently  the  U.  S.  Navy  is  using  the  XBT  system  manufactured 
by  the  Sippican  Corporation  of  Marion,  Massachusetts.  This  system  con¬ 
sists  of  an  expendable  shipboard  cannistcr  and  probe,  shipboard  launcher, 
and  a  strip  chart  recorder.  Data  so  obtained  an*  submitted  to  FNWC  and 
were  used  In  this  thesis. 
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The  original  Bureau  of  Ships'  specifications  (Danner ,  1966] 


for  an  XBT  system  are  given  in  Table  I. 


TABLE  I 


Temperature  range 
Temperature  accuracy 
Depth  range 
Depth  accuracy 

Ship's  speed  when  launching 


28F  to  96F 
+  0.4F 
0  to  1000  ft 

+  lit  or  15  ft,  whichever 
is  greater 

+  63^  of  a  temperature 
step  after  3  ft  depth 
change 


The  expendable  Sippican  probe  currently  used  has  a 
thermistor  sensor  in  its  nose.  The  probe  ia  connected  electrically 
to  the  recorder  by  two  conductor  wire 3  which  pay  out  from  both  the 
probe  and  the  canister.  The  can.iister  in  the  launcher  is  connected 
to  the  recorder  by  a  cable.  There  are  no  splices  in  the  wire.  The 
probe's  weight  is  a  critical  factor  in  determining  its  fall-rate. 

The  Sippican  system  was  found  to  be  within  +  3  gms  both  with  and 
without  its  Internal  copper  wire.  (The  deviation  was  fror  0.3 $  to 
0.64*  of  the  total  weight  of  the  probe.)  The  observed  average  fall- 
velocity  of  the  probe  was  20.3  +  0.2  ft/sec,  A  probe,  with  its 
Internal  wire,  had  an  average  terminal  velocity  of  19.9  +  0.7  ft/sec. 
Terminal  velocity  was  reached  at  18.0+0.5  ft  [Gouzie,  Sanders  and 
Littlehale ,  1966]  , 

The  present  XBT  strip  chart  recorder  is  a  Mark  2A  model 
using  pressure  sensitive  paper  for  recording  the  thermal  profile.  The 
paper  moves  veitlcaily  on  the  recorder  at  a  rate  of  3.25  in/min.  In  one 
minute,  the  XBT,  falling  at  20  ft/sec,  will  have  reached  a  depth  of  1200 
ft.  Thus,  1200  ft  will  be  recorded  on  3,25  inches  of  chart  paper,  or 
O.U027  of  an  inch  of  chart  paper  is  used  for  every  foot  of  fall. 
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The  chert  paper  records  horizontally  from  28P  to  96F  in 
7.0  inches,  or  about  0.103  of  an  inch  of  chart  paper  per  1.0F.  Table 
II  shows  the  resolution  capabilities  of  the  paper  using  different 
recording  Intervals. 


TABLE  II 

Resol ition 

Chart  Spacing  (Inches)  Temperature  (F)  Depth  (ft) 

0.01  0.0971  3.690 

0.005  0.0476  1.045 

0.001  0.0097  .369 

The  response  time  for  the  thermistor  is  110  milliseconds 
nominal  and  130  milliseconds  maximum  [Choate,  1970  and  Demeo,  1968]. 

A  time  constant-temperature  change  plot  shows  that  98.2  percent  of  a 
temperature  change  has  been  recorded  in  10.4  feet  of  fall,  83  percent 
recorded  in  5.2  feet  of  fall,  68  percent  in  two  feet  of  fall  and  35 
percent  in  one  foot  of  fall  (Figure  24). 

The  thermistor,  because  of  its  finite  sice,  has  thermal 
inertia  which  causes  it  to  act  as  «  low  pass  filter.  As  a  result, 
fine  structure  is  lost  or  at  best  smoothed  out.  The  original  trace 
could  be  reconstructed  if  the  filtering  characteristics  (of  both  phase 
and  amplitude)  of  the  thermistor  were  known.  These  could  be  obtained 
experimentally. 

b.  Data  Handling 

Present  procedures  for  handling  XBT  data  are  indicated  by 
the  log  sheets  which  accompany  the  XBT's  in  the  shipping  container. 
First,  a  BATHY  or  BAXBT  message  is  required  to  be  submitted  to  the 
appropriate  addresser  using  a  PRIORITY  precedence  message  [National 


41 


Arsunpt ions  : 


Q)  (IV)  »*ut<43  »in5BJ»dui»x 


42 


Figure  24.  Teaperature  Change  aa  Functions  of  Tine,  Tina  Constants,  and  Depth  for  Various 
Teaperature  Steps. 


Oceanographic  Data  Center,  1968,  1968a].  Upon  receipt  of  the  BATHY  or 
BAX8T  message  the  receiving  agency  uses  this  as  synoptic  data  for  sub¬ 
jective  ocean  forecasts  in  his  locality,;  The  BATHY  or  BAXBT  is  for¬ 
warded  through  the  Naval  Environmental  Data  Network  to  FNWC  in  Monterey. 
FNWC  uses  the  message  reports  to  update  its  ocean  fields  for  objective 
forecasting  in  various  oceanographic  and  meteorological  products. 

Upon  completion  of  a  cruise  the  shin  which  took  the  XBT 
cast  mails  the  original  trace  and  log  sheet  to  FNWC  for  digitizing. 

The  original  XBT  trace  is  digitized  for  use  by  FNWC,  and,  when 
requested,  the  results  of  the  digitizing  are  provided  to  the  XBT 
submitting  activity.  Upon  completion  of  digitizing,  the  original  XBT 
traces  are  forwarded  to  the  National  Oceanographic  Data  Center  (NODC; 
for  storage  in  their  archives.  NODC  also  receives  a  copy  of  the  digital 
data  on  tape.; 

Present  procedures  at  FNWC  use  DIGITBT  and  several  other 
specialized  XBT  data  handling  programs.  New  XBT  digitized  data  are 
stacked  on  the  XBT  Master  File  (presently  at  a  rate  of  about  500  traces 
per  week) . 

B.  OTHER  SENSORS 

Sensors  for  measuring  microthermal  structure,  other  than  those 
discussed  elsewhere  in  this  thesis,  are  numerous. 

From  a  fixed  platform,  like  the  NEL  Oceanographic  Tower,  thermistor 
strings  an  i  arrays  are  used.  Bottom  devices,  such  as  thermistors, 
attached  to  the  shore  or  a  tower  by  a  wire,  have  been  employed  to  give 
continuous  temperature  records. 
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Prom  submarines,  thermopiles,  thermocouples,  therm. stors,  platinum 
resistance  thermometers,  and  thin  film  thermometers  have  been  employed. 
Almost  any  oeeanograph ic  electronic  or  mechanical  thermometer  could  be 
used  on  a  submarine. 

C.  DIGITIZING  PROCEDURES 

At  the  present  time,  all  XBT  traces  are  submitted  to  Fleet  Numerical 
Weather  Central  (FNWC)  in  Monterey,  California,  the  only  digitising 
center  for  these  strip  chart  analogue  traces  (National  Oceanographic 
Data  Center,  1966,  1966a;  and  Bauer,  1969]. 

At  the  Point  Pinos  Annex  of  FNWC,  two  Calms  Model  480  flat-bed 
digital  recorders  are  used  to  digitise  XBT  traces  in  0.01  inch 
increments.  These  digital  recorders  are  operated  in  a  mode  which  makes 
the  recorded  intervals  dependent  on  the  tracing  speed.  The  Calma 
output  recorded  on  digital  magnetic  tape  containing  binary  mode 
characters  at  a  density  of  556  bits  per  inch  (BPI)  with  even  parity 
and  variable  length  records. 

FNWC  program  DIGITBT,  written  by  Lieutenant  N.  L.  Perkins,  USN, 
in  FORTRAN  IV  for  FNWC's  CDC  6500  computer,  Is  used  currently  for 
the  reduction  of  data  from  XBT  traces.  FNWC  tape  input,  tape  output, 
and  some  other  subroutines  of  DIGITBT  are  written  in  COMPASS  language. 

In  order  to  analyse  microthermal  structure,  DIGITBT  was  extensively 
modified.  The  new  program,  referred  to  as  M1CR0XBT,  is  included  as 
Appendix  A.  The  basic  prwc**urts  for  DIGITBT  and  MTCROXBT  programs 
are  the  same.. 

DIGITBT  an>i  MICR0X8T  programs  both  convert  the  X  and  Y  measure¬ 
ment  .  of  the  digit izei  into  depth  (X)  versus  temperature  (Y)  values. 


Corrections  are  made  for  non linear! ties  produced  by  the  XBT  recorder  in 
both  the  depth  end  temperature  directions.  All  shipboard  XBT  recorders 
are  assumed  to  produce  error-free  traces.  The  62. OF  temperature  line 
drawn  on  the  trace  when  the  probe  is  loaded  is  used  as  a  baseline  for 
digitizing.  No  correction  is  made  if  this  baseline  temperature  line 
is  not  the  same  as  the  62. OF  line  printed  on  the  trace  paper.  An 
eye-interpolated  value  of  the  XBT-drawn  baseline  temperature  on  the 
printed  temperature  grid  is  noted  as  a  part  of  the  identification 
record. 

The  rigid  procedures  for  digitizing  require  that  exactly  one 
identification  (ID)  record  and  two  digitized  trace  records  for  each 
XBT  be  recorded  on  the  tape.  The  ID  record  information  contains  the 
data  (day,  month,  year);  time  (hour  and  minute);  latitude  and  longi¬ 
tude  (degree,  minute,  hemisphere);  baseline  temperature  (to  nearest 
0.1C  or  0.1F  depending  on  the  grid  type);  printed-grid  type;  and  maxi¬ 
mum  depth  for  the  probe  (coded  as  15  for  1500  ft  XBT  or  25  for  che 
2500  ft  XBT).  Each  trace  record  contains  pairs  of  inrrements  (t$ C,  AY) 
read  from  the  baseline  to  the  true  surface  temperature,  a  point  "flag" 
and  then  pairs  of  values  (AX,  AY)  along  the  trace  to  a  cut-off  depth. 

Each  trace  is  digitized  twice.  The  complete  procedure  is  described 
by  Bauer  [1969], 

To  process  digitized  XBT  data  using  the  DIGITBT  program,  the  two 
trace  records  are  compared  point-by-point  in  a  deviation  test.  If 
the  X  and  Y  paired  values  (i.e.,  depth  and  temperature)  differ  by  more  than 
0.035  inches,  DIGITBT  rejects  the  whole  XBT  record  and  the  XBT  trace  must 
be  completely  re-digitized.  If  the  XBT  record  passes  the  deviation  test, 
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Temperature  unit*  arc  converted  from  digitised  inches  to  F  by  a  tempera* 
ture  table  and  from  7  to  C  by  formula.  These  depth  and  temperature 
tables  were  constructed  to  account  for  the  nonlinearities  between  the 
recorder  and  the  actual  trace.  Values  of  depth  are  rounded  to  the 
nearest  m  and  temperature  to  the  nearest  0.0 1C  in  D1GITBT.  MICROKBT 
has  eliminated  the  ft  to  m  and  P  to  C  conversions.  Depths  are  rounded 
to  the  nearest  0.1  ft  and  temperature  to  the  nearest  0.001F. 

After  the  conversion  from  ft  to  the  nearest  m,  DIGITBT  removes, 
for  a  second  time,  any  of  the  rounded-off  values  which  are  nonlinear 
between  two  neighboring  points.  If  depth  values  appear  as  duplicates, 
one  value  is  omitted  in  the  output  listing.  In  MICROXBT,  there  are  no 
duplicate  values  because  0.1  ft  in  depth  is  well  beyond  the  hundredth 
of  an  inch  digitising  capability  of  the  digitizer.  Thus,  by  reading 
depths  in  tenths  of  a  ft  and  temperature  in  thousandths  of  a  degree  F, 
MICROKBT  looses  no  values  from  duplicate  depths,  and  the  possibility  of 
introducing  false  linearities  is  reduced  considerably.  For  these  rea- 
sons,  MICROKBT  gives  a  much  finer  digitized  record  than  DIGITBT. 

The  DIGITBT  output  tape  is  written  in  binary  for  the  CDC  6500;  the 
MICROKBT  output  tape  is  written  for  the  IBM  360/67  at  the  Naval  Post¬ 
graduate  School  (NPS).  Appendix  0  contains  the  job  control  cards 
(JCL)  required  for  using  a  CDC  6500  output  tap'  from  MICROXBT  with 
external  BCD  at  a  density  of  556  BPI  as  in  input  to  a  separate  micro- 
thermal  analysis  program  at  NPS. 


IV  MODELS  FOR  EXAMINING  MICROTHERMAL  SUBCTOTI 


All  models  and  methods  used  in  this  thesis  required  the  output  tape 
from  MICROKBT  to  be  used  as  a  direct  input  on  the  seven  channel  tape 
unit.  The  XBT  ID  record  information  was  read  first  and  if  the  recorded 
XBT  was  from  a  desired  geographical  region  the  data  (depth  and  tempera- 
ture)  were  analysed  for  microthermal  structure  at  WPS. 

A  "best"  value  for  the  temperature  structure  was  found  by  using  the 
depth  inflection  points  from  MICROKBT.  For  this  "best"  array,  tempera¬ 
ture  was  linearity  interpolated  to  a  value  for  every  foot  of  depth.  The 
mixed  layer  depth  (MLD)  was  determined,  using  inflection  points  from 
MICROKBT,  as  that  depth  immediately  above  the  depth  where  the  temperature 
was  2. OF  different  from  that  at  the  surface  [Naval  Heather  Service,  1967]. 

Several  microstructure  models  were  tried  in  order  to  perform  an 
objective  analysis  of  the  oceanic  microthermal  structure.  The  use  of 
a  numerical  temperature  derivative  model  was  found  to  give  the  most 
satisfactory  presentation  of  the  microthermal  structure. 

A.  DERIVATIVE  MODEL 

1-  Mathematical  Model 

The  most  significant  feature  of  the  thermal  structure  is  the 
thermal  gradient  defined  as  the  derivative  of  the  temperature  with 
respect  to  depth  (dT/dZ).  By  plotting  the  gradient  versus  depth  (Z), 
the  thermal  structure  of  the  water  column  is  emphasised.  Furthermore, 
a  parameter  such  as  the  thermal  gradient  is  related  to  stability, 
entrophy  generation,  turbulent  heat  transport ,  and  the  Brunt -Vi lea la 


2.  laotrimental  Methods 


There  are  two  methods  for  experimentally  determining  tha  thermal 
gradient  presently  la  use. 

a.  Electronic 

A  resistance  thermometer  in  an  electronic  bridge  peasing 
through  water  produce*  a  signal  which  is  proportional  to  temperature  as 
a  function  of  tine.  By  differentiating  this  signal  with  respect  to  tine 
(which  is  analogous  to  differentiating  with  respect  to  distance),  the 
values  of  gradient  versus  distance  the  probe  travels  arc  fouul. 

This  method  was  successfully  employed  by  Grant,  Moll let  and 
Vogel  [1968],  in  the  horizontal  direction  for  temperature  and  turbulence 
measurements.  A  thin- film  resistance  thermometer,  mounted  on  the  bow  of 
a  submarine,  measured  temperature  versus  horizontal  distance  traveled. 
The  time  derivative  was  taken  to  produce  a  direct  plot  of  dT/dX,  the 
horizontal  thermal  gradients. 

Osborn  [1969]  used  a  thermistor,  capacitively  coupled  to 
two  differentiating  amplifiers,  to  measure  the  vertical  temperature 
gradient  with  his  free-falling  gradlometers .  His  resulting  output 
(Figure  13)  shows  a  gradient  versus  depth  record  for  the  North  Pacific 
off  the  San  Diego  Trough. 

b.  Two-probe  System 

Temperature  gradients  may  be  found  by  using  two  vertically 
spaced  thermistors.  The  temperature  difference  is  a  measure  of  the  mean 
gradient  between  them.  Such  a  gradient-meter  was  used  by  Woods  [1968, 
1968a]  who  obtained  gradient  at  various  depths  (Figures  7,  8,  and  9). 
This  method  is  limited  to  the  resolving  power  of  the  thermistors  and  is 
a  function  of  thermistor  spacing. 


3.  Bigtrlcil 


The  numerical  method  ia  bated  on  th«  darivatlva,  dT/dZ,  for  the 
gradient  but  usas  the  finite  difference  approximation.  The  darivatlva 
of  the  temperature  at  any  depth  T(Z^)  with  respect  to  depth  Z  is  defined 
as: 


A.  u.  Y<*i>  ’  Y<*> 
dt  Zj-Z  Zt  "  z 


<5> 


which  can  be  shewn  to  equal: 


nr  m,±m  :.m 

urt>  &  * 


(6) 


For  this  thesis i  At  •  1  foot,  which  yields  the  finite  difference  approxi¬ 
mation  for  the  vertical  temperature  derivative: 

(7) 


B.  MXfCILUMOUS  MMMOM 

1  MltiHl  liiUlUKiifltt 

For  this  nethod,  a  standard  or  snooth  curve  was  required. 

Several  standard  curves  were  Investigated.  The  variance  and  bias  were 
calculated  between  the  standard  and  the  "best"  curves.  The  standard  was 
found  by  computing  the  temperature  values  for  a  series  of  specified  depths. 

The  first  standard  curve  investigated  used  MLD  and  the  HIWC 
significant  levels  (surface,  100,  200,  300,  400,  600,  800,  and  1200  ft). 

The  temperature  values  were  found  using  the  "base"  temperature  values  at 
those  depths.  Hast ,  an  array  of  temperature  values  was  found  at  foot 
ifc  cements  t>»  linear  interpolation  between  the  temperature  values  at  the 
*  <§)•>!  fie  »»*  .*vss>i*  *t  wy*  yh » »  <***  Se«  the  FWbC  #t§#s?Je*«t 


l*v*l  curv*.  Calculations  of  bln  and  variance  did  not  ravaai  any  cor* 
relation  with  tba  microthermal  structure  easily  visible  on  the  X1T  trace. 

The  second  standard  curve  Investigated  used  the  temperature*  at 
the  standard  oceanographic  station  depths  (surface,  10,  ?0,  30,  90,  70, 
100,  150,  200,  250,  300,  400,  500,  snd  §00  «).  Since  these  depths  are 
nitric,  a  linear  interpolation  was  done  between  the  two  neighboring 
value*  on  the  "best"  curve.  Values  of  temperature  for  every  foot  were 
found  by  linear  interpolations  between  the  standard  dapth  temperature 
value*.  Again,  no  significant  correlation  w,  tad. 

Tha  third  standard  curva  usad  than  standard  depth  temperature 
valuta  but  uaed  thaaa  standard  dapth  valuta  as  an  input  to  tha  IflliffF 
subrouting  [Dtnnsr,  1969).  This  subroutine  uses  tha  four  standard  dapth 
points  nalghboring  tha  dapth  in  question  to  maha  an  "average"  curva  from 
an  upward  and  downward  parabolic  interpolation  and  a  linear  interpolation. 
This  schema  produces  a  smoothing  of  points  on  the  curve.  This  technique 
is  usually  usad  for  interpolating  standard  dapth  oceanographic  data 
[Hattray,  1965'.  These  values,  after  computing  tha  biaa  and  variance, 
again  gave  no  correlation  with  the  microthermal  structure. 

The  LGINTP  subroutine  wee  not  used  with  the  PNUC  significant 
level  values  for  smoothing  since  the  VNWC  fields  assume  linear  inter¬ 
polation. 

Other  standards  considered,  but  not  investigated,  included  a 
Gaussian  standard,  using  tha  technique  for  approximating  tha  temperature 
structure  proposed  by  Crosfils  [1968]  and  an  exponential  standard. 

2.  Gra4ieM^»tributi<^Mftthgd 

For  investigating  the  temperature  struct ure,  the  "heat"  tempera* 
tore  date  were  sepetste^  ini  *«#  PKIST  s,v*  i«  f  ?•<  h 


significant  level,  the  temperature  change  fro*  one  depth  to  the  next 
were  calculated  and  then  sorted  in  numerical  ascending  order  with  the 
IW  subroutine  SHSORT.  With  these  values  of  temperature  change  per  foot 
or  gradient  (P/ft),  the  temperature  structure  was  now  characterised  by 
the  individual  thermal  gradient  values  and  a  total  temperature  change  AT. 
Because  of  the  range  of  temperature  values  in  the  theraocline,  this  region 
of  the  XBT  trace  was  chosen  for  developing  this  method. 

The  capabilities  of  the  XBT  indicate  a  gradient  interval  of  0.02 
F/ft  is  easily  obtained  with  a  high  degree  of  confidence  [Arthur  D. 

Little,  1966].  Using  0.02  F/ft  for  the  class  interval,  the  frequency 
distribution  of  gradients  was  plotted  for  various  XBT  traces  by  signifi- 
cant  levels.  As  expected,  the  greatest  range  in  gradient  values  usually 
occurred  in  the  significant  level  just  below  the  mixed  layer  where  the 
range  of  temperature  values  AT  was  the  greatest. 

This  method  appears  to  give  a  good  indication  of  the  variability 
of  the  thermal  structure  in  the  depth  intervals  Investigated.  Several 
examples  of  the  different  types  of  gradient  distributions,  using  histo¬ 
grams  and  a  cumulative  frequency  diagram,  are  shown  in  Figures  25  through 
31.  The  actual  analysed  XBT  is  also  shown  on  the  left  of  each  figure. 
Looking  at  data  off  San  Diego,  the  first  three  gradient  distribution 
diagrams  (Figures  25,  26  and  27)  show  fairly  similar  cumulative  fre¬ 
quency  distributions.  The  next  two  (Figures  28  and  29)  show  nearly 
identical  distributions  for  water  inside  the  Gulf  Stream.  The  next  twc 
gradient  distributions  (Figures  30  and  31)  show  water  in  the  Gulf  Stream'.; 
frontal  region.  These  last  two  distributions  are  not  identical  for  many 
but  'rtaloSv  b*i«us#  the  positive  region  in  the  second 
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25.  Gradient  Distribution  for  the  100  to  200  Foot  Interval  for  a  Location  100  na  West  of 
San  Diego. 
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Figure  28.  Gradient  Distribution  for  the  100  to  200  Foot  Interval  for  Water  Inside  the  Gulf  Stream. 
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30.  Gradient  Distribution  for  the  1(0  to  200  Foot  Interval  for  Water 


100  to  200  Foot  Interval  for  Water  in  the  Gulf  Streaa' 


distribution  plot.  These  results  indicate  that  tht 
may  be  water  mass  related.  However,  this  could  not 
the  time  ailoted  for  this  thesis  but  la  recommended 
investigation . 


small  scale  structure 
be  fully  examined  in 
for  further 
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V.  APPLICATION  OF  NUMERICAL  DERIVATIVE  MODEL 

This  chapter  describes  the  application  of  the  numerical  derivative 
model  or  method  (NDM)  to  several  cases  of  real  and  assumed  data. 

A.  BASIC  CASES 

Graphs  are  used  to  show  the  features  of  the  NUN.  For  each  station, 
temperature  versus  depth  (T-Z)  and  temperature  gradient  versus  depth 
(G-Z)  plots  were  drawn  by  a  CalComp  plotter  using  the  computer  program 
with  the  DRAW  subroutine  in  Appendix  B.  The  scales  are  constant  for  all 
plots.  The  grid  squares  were  drawn  for  100-ft  depth  increments  from  the 
surface  to  a  depth  of  900  ft.  On  the  T-Z  plots  (reporduced  digitized  XBT 
data),  temperature  was  scaled  at  10P  per  grid  line.  On  the  G-Z  plots,  the 
gradient  was  scaled  to  O.lF/ft  per  grid  Interval,  with  chi  axis  offset  so 
that  two-thirds  of  the  plot  was  allocated  to  negative  gradients.  This 
scaling  was  chosen  so  that  each  hundredth  of  an  inch,  the  resolving 
power  of  the  CalComp  plotter,  represented  one  foot  of  depth  and  0.01F 
of  temoerature  or  O.OOlF/ft  of  temperature  gradient.  The  plots  included 
here  have  been  reduced  in  size,  but  the  relative  proportions  remain  the 
same . 

1.  Linear  Gradients 

For  a  temperature  distribution  with  linear  thermal  gradients, 
the  G-Z  plot  is  characterized  by  steps  indicating  a  constant  gradient 
during  a  temperature  change.  For  this  thesis  the  strength  of  a  gradient 
is  defined  as  the  absolute  value  of  the  gradient  (in  F/ft  units)  and  is 
described  as  strong  when  greater  than  0.3F/ft  or  weak  when  less  than 
0.05F/ft.  The  term  signed  strength  is  defined  as  the  true  gradient 
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value  and  will  thus  distinguish  absolute  and  true  gradient  values.  The 
magnitude  of  a  gradient  is  defined  as  the  vertical  extent  of  the  gradient 
measured  in  feet.  The  product  of  magnitude  and  signed  strength  gives  the 
total  temperature  change  for  the  extent  of  the  gradient. 

a.  Layer  Structure 

Figure  32  shows  a  two- layer  temperature  structure  model  with 
a  T-Z  plot  on  the  left  and  a  G-Z  plot  on  the  right.  The  strength  of  the 
upper  layer  is  0.05F/ft  with  a  magnitude  of  200  ft.  The  lower  layer  is 
weaker  (0.02F/ft  compared  to  0.05F/ft)  than  the  upper  step  but  the  magni¬ 
tude  is  greater  (600  ft  to  200  ft). 

For  the  multi-layered  model  (Figure  33)  several  linear 
gradient  steps  are  shown.  In  the  gradient  region  between  70  and  90  ft, 
the  gradient  is  strong  and  goes  off  scale  (Figure  33b). 

b.  Step  Structure 

A  step  structure  is  characterized  by  thin  strong  gradient 
regions  The  model  (Figure  34)  shows  clearly  that  the  four  steps  bet¬ 
ween  60  and  100  ft  are  stronger  than  any  other  steps  shown.  These  four 
steps  also  have  the  smallest  magnitude  (five  ft)  in  the  figure.  The 
steps  at  200,  300,  and  400  ft  all  have  the  same  strength  but  increase 
in  magnitude  with  depth. 

The  features  of  the  T-Z  trace  (Figure  34a)  are  easily 
described  by  th<»  G-Z  model.  The  product  of  signed  strength  and  magnitude 
gives  che  total  tempera «.«..re  change  from  the  surface  temperature  to  a 
given  depth  for  a  step-structure  temperature  profile.  The  temperature 
at  any  depth  could  then  be  found  if  the  sea  surface  temperature  was 
known.  Thus,  with  magnitude,  signed  strength  and  sea  surface  temperature, 
any  layered  or  stepped  thermal  profile  could  hr  reconstructed. 
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Flgurt  33.  Multllayar  Modal 


Gradient  (G) 


2.  Continuous  Curving  Gradients 


Since  t'.nite  values  are  used  for  both  depth  and  temperature,  a 
truly  continuous  curve  cannot  be  drawn.  The  spacing  by  the  CalComp 
plotter  of  every-  hundredth  of  an  inch  can  cause  variations  of  this 
size.  Computed  values  of  depth,  temperature  and  temperature  gradient 
were  provided  as  inputs  to  a  CalComp  plotter  at  every  hundredth  of  an 
inch  to  create  a  quas i -continuous  curve,  at  least  to  a  resolution  of  the 
plotter.  It  is  this  quasi -continuous  field  which  accounts  for  the  0.01 
inch  steps  encountered  in  curves  drawn  by  the  plotter  for  any  truly  con¬ 
tinuous  curve.  This  effect  is  also  present  whenever  the  plotter  is 
drawing  any  line  not  directly  in  the  X  or  Y  axis  directions.  Accordingly, 
any  steps  shown  of  this  size  could  be  considered  as  possible  adjustments 
made  by  the  plotter  or  data  to  adjust  to  the  hundredth  of  an  inch  grid. 

Two  examples  of  continuous  gradients  are  given  to  show  the 
resulting  distributions  of  gradient  with  respect  to  depth, 
a.  Gaussian 

An  example  of  a  Gaussian  profile  (Figure  35)  shows  a  normal 
distribution  of  temperature  with  depth  for  the  function,  like 


T(Z)  -  80 


20  r 

ft  j 


-i«2 

2 

e  dx 


(8) 


where:  x  »  a  dummy  variable 

80  ■  sea  surface  temperature. 

The  G-Z  curve  highlights  the  changes  in  the  T-Z  curve.  The  increase  in 
gradient  strength  is  very  clear,  reaching  a  maximum  200  feet,  whereas 
the  inflection  point  of  the  T-Z  plot  is  difficult  to  find  by  visual 
inspection.  The  summation  of  signed  strength  values  at  every  foot  of 
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depth  from  the  surface  down  gives  a  total  temperature  change  which  is 
analogous  tc  the  product  of  magnitude  and  signed  strength  with  all 
magnitudes  set  equal  to  one.  With  this  summation  and  a  given  sea  surface 
temperature  (SST) ,  a  continuous  temperature  profile  can  be  reconstructed, 
b.  Exponential 

Using  an  exponential  expression  for  the  temperature  distri- 


but ion 

1(2)  -  60  -  20^/m0 


(9) 


where  : 


60  =»  SST, 

the  T-Z  (Figure  36a)  and  G-Z  curves  (Figure  36b)  appear  as  compliments 


of  each  other.  The  derivative  of  T(Z)  with  respect  to  Z  gives: 

-nZ/ 1000 


dT(Z) 

dZ 


TTe 

50 


(10) 


which  is  the  equation  of  the  G-Z  plot. 

3  Superpos i t ions 

In  the  example  of  layers  superimposed  on  an  exponential  model 
(Figure  37),  the  conditions  for  the  exponential  (9)  hold  for  the  depth 
intervals:  surface  to  40  ft,  220  to  340  ft  and  460  to  800  ft,  with 
various  layers  at  the  remaining  depths.  Isothermal  layers  exist  below 
the  weak  surface  expenencial  gradient  layer  and  also  in  the  140  to  180  ft 
depth  layer.  From  the  T-Z  plot,  the  layers  at  340  and  ^o0  ft  are  not  as 
evident  as  they  are  in  the  G-Z  plot,. 

4 .  Invers ions 

A  model  with  two  inversions  and  several  linear  steps  (Figure  38) 
shows  the  positive  gradient  regions  clearly.  The  second  deeper  inversion 
n  clearly  a  weak  gradient,  but  it  is  still  discernible  in  Figure  38b. 
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Figure  37.  Exponential  Model  with  Layers  Superimposed. 
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With  inversions,  summing  the  signed  strength  values  for  each  foot  gives 
the  temperature  change  from  the  surface  and  adding  this  value  to  SST 
yields  the  actual  temperature  at  any  depth. 


B.  OBSERVATIONAL  CASES 


Several  different  oceanic  regions  were  selected  to  establish  fea¬ 
tures  of  this  model.  The  oceanic  regions  explored  were  along  the  Atlantic 
Coast^in  the  Andaman  Sea,  Eastern  North  Pacific,  and  the  Gulf  Stream. 

1.  Atlantic  Coastal  Region 

Stations  in  Massachusetts  Bay  occupied  during  December  1969 
were  used.  The  stations  were  within  a  radius  of  13  miles.  Water  depths 
were  between  280  and  420  feet.  Lack  of  strong  gradients  in  all  seven 
traces  shows  that  the  water  was  probably  well  mixed  (Figures  39  through 
42).  A  positively  signed  strength  in  the  upper  300  feet  was  common  in 
the  traces.  Inversions  (Figure  40)  are  quite  clear  in  both  the  T-Z  and 
G-Z  plots.  A  portion  of  the  actual  XBT  trace  is  shown  for  three  of  the 
stations.  In  these  three  cases,  there  is  excellent  visible  correlation 
between  the  actual  XBT  trace  and  the  reconstituted  T-Z  plots  which  shows 
that  MICROXBT  can  resolve  very  small  temperature  features.  In  no  case 
was  there  an  isothermal  layer  greater  than  140  feet  thick. 

2 .  Andaman  Sea 

In  the  Andaman  Sea,  near  the  Straits  of  Maalacca,  a  series  of 
15  XBT  stations  from  the  18  to  20  July  1969  period  were  used.  The 
stations  were  within  a  radius  of  12  miles.  Water  depths  varied  between 
290  and  380  feet.  This  oceanic  region  could  be  considered  somewhat 
similar  to  the  Atlantic  Coastal  region  above. 
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Figure  39.  T-Z  and  G-Z  Plots  for  Atlantic  Coastal  Region  6  and  7  December  1969. 
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Figure  40.  T-Z  and  G-Z  Plots  for  Atlantic  Coastal  Region  7  December  1969. 
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Figure  41.  T-Z  and  G'z  Plots  for  Atlantic  Coastal  Region  9  and  10  December  1969 
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Pigure  42.  T-Z  and  G-Z  Plots  for  Atlantic  Coastal  Region  11  December  1969. 


The  T-Z  plots  (Figures  43  through  50)  contain  a  mixed  layer  to 
200  or  more  feet  except  for  one  observation  at  1200Z  on  the  19th  (Figure 
48b).  The  mixed  layer  transients  are  marked  with  spikes  in  the  G-Z  plots. 
Figure  50  shows  transients  with  the  greatest  magnitude  in  the  upper  100 
feet.  The  seasonal  thermocline  is  characterized  by  two  steps  (Figure 
44,  for  example).  The  G-Z  plots  in  most  cases  show  this  second  step, 
averaging  at  a  depth  about  300  feet,  with  good  resolution.  Some  strong 
gradients  are  found  in  the  200  to  300  feet  depth  region  in  every  G-Z 
plot  except  in  Figure  48b.  This  figure  appears  to  indicate  flow  over  a 
hump  or  have  erroneous  depth  values.  The  fall-rare  of  the  XBT  probe, 
appearing  erroneously  in  the  T-Z  plot,  indicates  a  possible  malfunction  — 
of  the  probe  or  the  recorder  for  this  one  cast. 

3.  Pacific  Oceanic  Region 

A  series  of  XBT  traces,  taken  over  a  week  period  in  October  1968, 
in  the  deep  waters  (over  100  nm  off  the  West  Coast  near  San  Diego, 

California)  were  used  for  this  observation.  A  composite  T-Z  piot  of 
the  actual  XBT  temperature -depth  traces  is  shown  m  Figure  51.  This 
regio.i  is  sub-tropical  oceanic  (Tully,  19641.  Plots  of  T-Z  and  G-Z 

(Figures  52  through  65)  are  presented  from  the  surface  to  900  feet.  The - 

seasonal  thermocline  is  evident  in  all  plots.  The  strong  thermocline  is 
clearly  shown.  Some  transients  are  notable.  Few  inversion  layers  were 
present,  and,  generally,  they  were  found  only  in  the  near-surface  layer. 

The  G-Z  plot  is  characterized  by  strong  gradients  in  the  seasonal  thermo¬ 
cline  and  by  weak  gradients  in  the  main  thermocline,  each  separate  layer 
separated  by  isothermal  waters  from  one  to  forty  feet  thick. 

The  mixed  layer  contained  numerous  temperature  fluctuations  in 
most  cases.  Some  inversions,  in  particular  the  one  at  0600Z  on  14  October 


77 


Z(xl02)ft 

180769  1800  0600N  09815E 

G(xK)*1)  P/ft  T(xlOl)F 

•30*;  -003  -032  -COi  OjO  Qjj 


180769  2000  0600N  098156 


Figure  44.,  T-Z  and  G-Z  Plots  for  Andaman  Sea  Region  for  18  July  1969. 
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Figure  45.  T-Z  and  G~Z  plots  for  Andaman  Sea  Region  for  18  and  19  July  1969 
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Figure  46.  T-Z  and  G-Z  Plots  for  Andaman  Sea  Region  for  i9  July  1969. 
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Figure  47.  T-Z  and  G-Z  plots  for  Andaman  Sea  Region  for  19  July  1969. 
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Figure  48.  T-Z  and  G-Z  Plots  for  Andaman  Sea  Region  for  19  July  1969. 
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Figure  49.  T-A  and  G-Z  Plot*  for  Andaman  Sea  Region  for  19  July  1969. 
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Figure  50.  T-Z  and  G-Z  Plots  for  Andaman  Jea  Region  for  20  July  1969. 
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Figure  51.  Composite  T-Z  Plot  from  Actual  XBT  Traces  Down  to  900  feet 
for  Region  100  nm  off  San  Diego.  Data  are  given  in  time 
sequence  from  10  through  16  October  1968. 


101068  180Gz(xl°  )ft3246N  12043W 

Figure  52.  T-Z  and  GtZ  Plots  for  Eastern  North  Pacific  for  1800Z  10  October  1968. 
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Figure  53.  T-Z  and  G-Z  Plots  for  Eastern  North  Pacific  for  OftOOZ  11  October  1968. 
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Figure  54.  T-Z  and  G-Z  Plots  for  Eastern  North  Pacific  for  0000Z  12  October  1968. 
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Figure  55.  T-Z  ar.d  G-Z  Plots  for  Eastern  North  Pacific  for  0600Z  12  October  1968 


Figure  56.  T-Z  and  G-Z  Plots  for  Eastern  North  Pacific  for  1800Z  12  October  1968. 
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Figure  57.  T-Z  and  G-Z  plots  for  Eastern  North  Pacific  for  0000Z  13  October  1968. 
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Figure  58.  T-Z  end  G-Z  Plots  for  Eastern  North  Pacific  for  0600Z  13  October  1968. 
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Figure  59«  T-Z  and  G-Z  Plot*  for  Eaatirn  North  Pacific  for  1200Z  13  October  1968. 
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Figure  60.  T-Z  and  G-Z  Plots  for  Eastern  North  Pacific  for  1800Z  13  October  1968. 
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Figure  61.  T-Z  end  G-Z  Plots  for  Eastern  North  Pacific  for  0600Z  14  October  1968. 
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Figure  62.  T-Z  and  G-Z  Plota  for  Eastern  North  Pacific  for  1200Z  14  October  1968. 
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Figure  63.  T-Z  and  G-Z  Plots  for  Eastern  North  Pacific  for  1800Z  14  October  1968. 
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1968  at  120  feet  was  clear  on  both  the  ba thy thermogram  and  the  T-Z  plot. 
The  steps  are  quite  irregular  on  the  T-Z  plot  but  not  as  irregular  as  on 
the  G-Z  plot. 

4  Gulf  Stream  Region 

Thirteen  stations,  during  the  period  11-13  August  1969,  were 
analysed  for  the  Gulf  Stream  region  (Figures  66  through  80).  The  ship's 
track  crossed  the  Gulf  Stream  boundary  several  times. 

a.  Interior  of  Gulf  Stream 

For  the  interior  of  the  Gulf  Stream  (Figures  70  and  71), 
the  mixed  layer  is  distinct  and  a  typical  seasonal  and  main  thermocline 
are  shown.  A  "boot"  in  the  100  foot  layer  is  clearly  shown  by  positive 
gradient  values  in  several  of  the  figures  and  the  actual  XBT  traces. 

b.  Frontal  Region 

Strong  gradients  are  the  rule  for  this  region  and  they  are 
quite  distinct  in  both  the  T-Z  and  G-Z  plots.  The  contrast  between 
the  frontal  and  interior  regions  of  the  Gulf  Stream  are  quite  apparent. 
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Figure  66  T-Z  and  G-Z  Profile  for  Transects  across  Northern  Boundarv  of  the  Gulf  Stream.  Actual  XBT  T-Z  Trace  Inset 
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Pigurt  68,:  T-Z  if.d  G«Z  Profile  for  Trarnacts  acroi*  Northern  BoirMar\  'f  Mu*  Gulf  Stream,  Actual  X8T  T-Z 
Trace  Inaet  . 
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rigurt  69,  T-Z  and  G-Z  froflla  for  Transacts  across  Northarn  Boundary  of  tha  Gulf  Stream.  Actual  XBT  T-Z 
Tract  Insat, 


Fi|ur*  70.  T-l  «nd  O-l  frofiU  for  Tranaacta  actum  Northern  Boundary  of  the  Gulf  St rate.  Actual  XIT  T-Z 
traca  Inset. 
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Flgura  71.  T-Z  and  G*Z  Froflla  for  Tranaacta  acroaa  Xortharn  Boundary  of  tha  Gulf  St  roam..  Actual  XBT  T*Z 
Traca  Xnaat. 
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Figuro  74,  T*X  and  G-2  SProfllo  for  Transacts  across  Nor  thorn.  Boundary  of  tho  Gulf  Straam. 
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Figure  75.  T-Z  and  G-l  Frofile  for  Trantacta  across  Northern 
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of  the  Gulf  Stream. 


Ftjur#  76.  T-2  and  C*X  ftoflU  for  TranitUa  acroaa  Hortharn  Boundary  of  tha  Gulf  Straam. 
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rigut*  77.  T-Z  and  0*1  frofll#  for  Tranocti  icrou  Northtrn  Boundary  of  the  Gulf  Strom. 
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Pigur#  78*r  M  and  G-Z  Pro ft  I#  for  Transect*  acroas  Northern  Soundary  of  tha  Gulf  Stream,, 
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VI.  CONCLUSIONS 


The  results  give  a  measure  of  the  microthermal  structure  obtained 
with  an  XBT  by  the  NDM. 

A.  PHYSICAL  SIGNIFICANCE 

Microthermal  structure  statistics  have  not  been  adequately  explored. 
This  thesis  has  proposed  an  easy  method  to  observe  the  fine-s:alt;  thermal 
structure.  Further,  this  proposed  objective  analysis  scheme  does  not  use 
the  T-Z  values  but  rather  the  G-Z  values.  It  Is  the  G-Z  distribution 
vhich  is  important.  The  T-Z  ran  always  be  reconstructed  if  one  additional 
parameter,  the  SST,  is  known. 

1.  Internal  Waves  -  Vaisala  Frequency 

The  Brunt-Vaisala  frequency  is  the  maximum  frequency  of  oscil¬ 
lation  of  free  internal  waves.  It  may  be  calculated  from  the  density 
gradient,  or  from  the  parameters  which  define  density  gradient.  In 
isohaline  water,  the  NDM  gives  a  direct  measure  of  density  (or  stability). 
In  such  waters  it  could  be  used  to  set  an  upper  limit  on  the  maximum 
frequency  of  oscillation  of  internal  waves.  Put  I  her,  a  time  series 
study  would  clearly  illustrate  the  amplitudes  likely  to  be  encountered 
since  the  sheets  indicated  on  the  G-Z  diagram  are  readily  tdenclfied. 

A  correlation  study  of  G-Z,  Brunt-Vaisala  versus  depth  diagrams  would 
be  useful. 

2.  Oceanic  Fronts 

This  difficult  region  to  analyze  microscopically  could  possible 
be  better  explored  by  using  the  NDM  for  the  microthermal  processes.  The 
NDM  explores  the  continuously  changing  conditions  in  frontal  regions. 


117 


Frontal  regions  are  usually  Indicated  by  extremely  complex  G-Z  plots 
which  illustrate  the  extreme  gradients,  both  positive  and  negative, 
encountered.  Based  on  the  admittedly  small  number  of  cases  examined, 
there  appears  to  be  a  tendency  of  the  layers  to  maintain  their  character¬ 
istics  and  slide  by,  over  or  under  each  other  rather  than  mix. 

3 .  Sound  Propagation 

The  effect  of  mlcrothermal  structure  on  sound  propagation 
require  hatter  statistical  knowledge  of  the  oceans.  The  NDM  gives 
emphasis  to  the  strongest  gradient  regions. 

The  mlcrothermal  effects  on  index  or  refraction  variations 
(density  fluctuations)  must  be  determined.  Using  the  NDM  with  salinity 
variations  (from  an  STD  cast)  and  the  NDM  with  sound  velocity  variations 
(from  an  SVTP  cast)  taken  from  the  two  instruments  simultaneously,  cross- 
correlation  techniques  could  shed  some  light  on  the  unresolved  micro¬ 
processes  in  the  ocean.  In  particular,  what  are  the  exact  physical  and 
theoretical  correlations  among  the  salinity,  temperature,  and  sound 
velocity. 


B .  SUGGESTS?)  STUDIES 


Possible  studies  for  the 
developing  this  thesis. 

1 .  Mathematical  Models 
Various  mathematical 
because  each  method  had  some 
their  utility. 

a.  Analysis  Methods 
The  mlcrothermal 
as  a  "roughness"  factor,  the 


future  have  suggested  themselves  while 


models  attempted  could  be  explored  further 
merit.  Suggested  models  are  limited  by 


structure  recorded  by  XBTs  can  be  looked  at 
noise  level  of  the  water  structure.  Further 
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work  is  required  to  determine  such  a  factor  for  practical  operational  use. 
This  roughness  factor  should  be  computed  as  the  XBT  trace  is  being 
analyzed  at  FNWC  by  a  program  like  DIGITBT .  This  factor  could  be  one 
number  to  describe  objectively  the  microthermal  structure  variations  on 
the  G-Z  trace. 

Variations  in  the  NDM  should  be  used  with  the  variables  that 
have  been  arbit  srily  set.  For  example,  the  depth  interval  of  AT/ft 
at  one  foot  should  be  varied  t r»  find  the  effects  a  larger  or  smaller 
depth  interval  would  have  on  the  G-Z  profiles.  As  the  steps  become 
smaller,  the  exact  derivative  is  approached.  Larger  steps  should  smooth 
the  spikes  shown  in  the  G-Z  figures.  In  some  cases,  this  smoothed  G-Z 
profile  might  be  described  with  curve  fitting  techniques  as  suggested  by 
Boston  [1966]  and  Grosfils  [1968]. 

MICROXBT’s  digitizing  interval  of  0.01  inches  should  be 
changed  to  0.005  inches  and  the  program  modified  to  use  this  interval. 

This  would  yield  mure  detailed  thermal  structure  than  MICROXBT.  Obviously, 
one  cannot:  make  the  interval  finer  than  the  noise  level  of  the  sensor  or 
recording  instrument . 

b.  Physical  Processes 

The  G-Z  profiles  from  the  NDM  emphasize  the  changing  physical 
features  which  are  obscure  in  the  T-Z  profiler.  By  using  G-Z  profiles 
for  parameters  other  than  temperature,  clarity  in  determining  distri¬ 
butions  of  chemical  parameters  would  be  similarly  emphas  J.zed .  Limitations 
in  use  of  the  NDM  on  other  chemical  properties  should  be  further 
investigated. 

Water  mass  identification  is  possible  with  G-Z  profiles. 
Additional  work  is  required  in  this  area. 


119 


2.  Experiments  l  Work 


A  mlcrothermal  sheet  probe  could  be  designed  for  investigating  the 
chemical  end  physical  processes  inside  the  sheet  end  inside  the  layers. 

A  comparison  of  actual  physical  property  differences  should  be  made. 

A  modified  expendable  temperature  gradient  (XTG)  meter  could  be 
designed  using  two  thermistors  to  obtain  the  temperature  differences  in 
the  probe.  The  standard  XBT  system  equipeient  already  Installed  could  be 
changed  to  give  analog  G-Z  profiles  to  the  operational  user  without  any 
delay. 

e.  Time  Series 

The  XTG  meter  could  be  used  to  construct  G-Z  variability  at 
one  location  by  conducting  a  time  series  analysis,  nils  needs  to  be 
done  in  order  to  obtain  e  roughness  factor ,  an  objective  microstructure 
measuring  parameter. 

b.  Transect 

An  XTG  meter  could  be  used  in  transects  across  various  water 
messes  to  obtain  the  definite  characteristics  for  the  water  masses.  The 
XTG  could  be  used  to  follow  interesting  features  on  the  G-Z  plots.  More 
important,  however,  is  that  the  special  characteristics  of  microtherma 1 
structure  could  easily  be  determined. 

3.  Handling  of  Data 

When  the  microthermal  structure  roughness  factor  is  determined, 
this  factor  should  be  made  an  integral  part  of  the  records  for  that  XBT 
trace,  nils  Information  should  be  included  as  pert  of  the  ID  data. 

Aboard  Ravy  ships  the  AH/UNQ-7  tape  recorder  could  be  used  to 
obtain  magnetic  tape  records  of  the  ere  log  X1T  or  XTG  signals.  These 
tsp*»  mild  b#  t.ifeni t f with  the  tree*#  ftWC  for  further  sn»l> 


4.  Maw!  Application* 


lasults  from  nirrothermal  studies  show  that  high  frequency  sonar 
performance  as  influenced  by  the  microstructure  receive  the  greatest 
nic restructure  effects.  More  knowledge  of  the  microstructure  is  needed 
before  the  torpedo  and  mine-hunting  sonar  performance  could  be  improved. 
With  a  measure  of  the  microstructure,  the  performance  of  a  sonar  system 
could  be  predetermined. 

Tens  of  thousands  of  XBT  traces  are  on  file  at  FNWC.  Hew  data 
are  arriving  daily.  This  thesis  is  an  attempt  to  show  how  this  data  can 
be  used  both  operationally  and  scientifically.  Obviously  it  represents 
only  a  beginning  in  the  way  of  microstructure  analysis  and  operational 


utilisation. 


APPENDIX  A 


MICROKBT  PORTRAN  IV  PROGRAM  FOR  COC  6500 


PROGRAM  MICROXBT  CINPUT, OUTPUT, TAPE21 
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//(Job  name) 

J0B  (Acct.  no..  Job  no,.  Section),  ’(Your  name)' 

Note  1 

//  EXEC 

F0RTCLGP,REGION«1OOK,  TIME .  G0-2 

Note  2 

//F0RT.SYSIN 

DD  * 

FORTRAN  Source  Deck 

/* 

Note  3 

//G0.FTO6FOO1 

DD  SYS0UT-A,SPACE(CYL,6) 

Note  4 

//G0.FTO4FOO1 

DD  UNIT-2400- 1,V0L-SER-NPSXXX,  LABE  W,NL) , 

Note  5 

// 

DCB-(DEN-1,RECFM-F,LRECI^136,BLKSIZE-136,TRTCH-ET),  ”  5 

// 

DISP-(0ID,KEEP) 

Note  5 

//G0.FTO4FOO2 

DD  AFT-FT04F001 , V0L-SER-NPSXXX , LABEL(2 ,NL) , 

Note  6 

// 

DCB-(DEN-1,RECFM-F,IHECL-136,BLKSIZE-136,TRTCH-ET),  "  6 

// 

DISP-(0LD,KEEP) 

Note  6 

//C#. SYS IN 

DD  *  '  i 

Note  7 

DATA  Cards,  if  any 

/*  (Orange  Card) 

Notes 

1.  Standard  J0B,  green  card,  for  program  submission. 


2.  If  not  plotting  or  card  punching,  the  F0RTCLGP  can  be  changed 
to  F0RTCLG. 

3.  Must  be  a  white  card. 

4.  Optional  card  for  additional  storage  for  output. 

5.  The  XXX  following  NPS  is  the  NPS  tape  number.  These  three  cards 
are  for  reading  the  first  file  on  the  input  tape.  These  cards 
assign  input  unit  4  for  the  source  programs  READ  statements. 
Further  note  that  the  record  length  is  136  characters,  all  of 
which  may  be  read.  This  is  the  standard  record  length  from  the 
CDC  3600. 

6.  For  a  third  file  to  be  read,  these  three  cards  must  be  repeated 
with  F002  changed  to  F003  and  IABEL(2 ,NL)  changed  to  LABEL(3 ,NL) 
and  so  forth  for  additional  files.  Also  the  UNIT-2400-1  should 
be  changed  to  AFT-FT04F001  for  all  files  after  the  first  file. 

If  only  one  file  is  to  be  read,  then  these  three  cards  are  not 
required . 

7.  Not  required  if  there  are  no  input  data  cards. 
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JOB  CONTROL  LANGUAGE  (JCL)  FOR  DUMPING  FNWC  TAPE  DATA 
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//(Job  name) 
//DUMPTAPE 
//SYS PRINT 
/ /SYSUT1 
// 

//SYSUT2 
//SYS IN 


J0B  (Acct.  no.,  Job  no.,  Section),  ’(Your  name)'  Note  1 

EXEC  PGM-IEBPTPCH 
DD  SYS0UT=A 

DD  UNIT=24OO-1,DI3B=0LD,IABED»(,NL),V0LUME=SER=NPSXX,  Note  2 

DCB»(DEN=l,RECFM=F,LRECIJ«136,BLKSIZE»136,TR,rCH»ET)  Nofe  2 

DD  SYS0UT=A 
DD  *,DCB»BLKSJ,ZE=80 

PRINT  MAXFIDS=1,ST0PAFT“IOOO  Note  3 

RECORD  FIELDS=(120)  ‘  Note  4 


/*  (Orange  Card) 


Notes 


1.  Standard  J0B,  green  cavd,  for  program  submission. 

2.  If  more  than  one  file  is  on  che  tape  then  separate  two  cards  for 
each  file  identified  by  file  number  by  the  file  number  in  the  IABEL 
(X,NL)  where  S  is  the  file  numbur,  which  is  implied  1  when  omitted. 

3.  The  ST0PAFT  limits  the  number  of  lines  of  printed  output,  in  this 
example  1000  lines  is  used  before  printing  is  terminated. 

4.  The  field  is  recorded  at  the  maximum  amount  that  this  program  will 
dump  which  is  120  characters.  If  there  are  any  characters  between 
120  and  133  in  the  records  then  this  routine  will  not  dump  those 
characters . 
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